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NUMERICAL  STUDY  OF  PHASE  CONJUGATION 
IN  STIMULATED  BACKSCATTER  WITH 
PUMP  DEPLETION 


I.  INTRODUCTION 

Optical  wavefront  reversal  or  phase  conjugation  in  stimulated  backscatter  has  been  the  subject  of 
extensive  experimental1'15  and  theoretical14'22  investigation  over  the  last  ten  years.  Most  of  the  experi¬ 
ments  are  either  of  the  aberration  correction  (ABC)  type,1'9  in  which  the  conjugate  backscatter  reverses 
phase  distortion  induced  on  the  incident  beam  by  an  aberrator,  or  the  image  replication  (IMR)  type,10'13 
in  which  the  backscatter  retraces  an  image  formed  by  a  mask  (with  or  without  an  aberrator)  in  the 
incident  beam.  The  ABC  experiments  normally  record  only  the  far  field  of  the  aberrator,  whereas  the 
IMR  experiments  observe  at  a  point  conjugate  to  the  mask.  Although  the  existing  analytic  theories14'20 
can  explain  how  a  phase  conjugate  wave  is  generated  in  stimulated  backscatter,  they  have  made  no 
attempt  to  actually  model  these  experimental  configurations;  nor  have  they  realistically  included  pump 
depletion,  which  becomes  important  when  the  backscatter  is  significantly  above  threshold.17-913  (One 
analytic  treatment  has  included  pump  depletion,21  but  its  assumptions  and  limitations  are  not  applicable 
to  the  conditions  studied  here.) 

In  an  earlier  paper,22  image  replication  in  stimulated  Briilouin  scattering  (SBS)  of  coherent  beams 
was  studied  numerically,  using  a  steady  state  2D  propagation  code  (BOUNCE).  The  present  work  treats 
phase  conjugation  of  a  focused  aberrated  beam,  using  a  modified  version  of  BOUNCE  that  has  been 
extended  to  include  pump  depletion.  In  all  of  the  cases  studied,  pump  depletion  significantly  enhanced 
the  fidelity  of  the  phase  conjugation  process  by  inhibiting  spatial  gain  narrowing.  The  simulations  show 
that  the  far  field  distribution  of  light  scattered  back  through  the  phase  aberrator  exhibits  a  prominent 
axial  spike  closely  matching  that  of  the  incident  beam,  in  agreement  with  the  ABC  experiments.  How¬ 
ever,  the  near  field  intensity  exhibits  large  and  rapid  spatial  inhomogeneities,  even  when  the  conjuga¬ 
tion  fidelity  approaches  90%.  In  spite  of  such  structure,  the  backscatter  was  able  to  reproduce  a  rough 
image  of  large  scale  intensity  modulation  imposed  upon  the  incident  beam;  i.e.,  the  results  can  simulate 
even  those  IMR  experiments  in  which  phase  aberration  was  also  present.10- 11 

II.  THEORY  AND  CODE  DESCRIPTION 

Under  steady  state  conditions,  the  complex  pump  and  SBS  backscatter  amplitudes  £t0c,r)  and 
Es(x,z)  satisfy  the  parabolic  equations 

El  -  jg\Es\2EL  (la) 

Es-  y g \ElVEs  (lb) 

in  a  two  dimensional  cartesian  geometry.  Here,  k  -  2nn/\  is  the  magnitude  of  the  propagation  vectors 
(assumed  equal  for  the  two  waves),  x  is  the  transverse  coordinate,  and  g  is  the  coupling  coefficient  of 
the  Briilouin  medium,  which  is  contained  within  a  region  zt  ^  z  <  z2.  BOUNCE  solves  Eqs.  (I), 
assuming  the  aberrated  pump  wave  is  incident  at  z2,  while  the  backscatter  grows  from  a  small  counter- 
propagating  noise  wave  £s(x,ri)  introduced  at  z\.  (E.g.,  see  Fig.  1.)  Both  waves  are  assumed  to  vanish 
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at  (and  near)  the  transverse  boundaries  ±X/2.  The  calculations  were  carried  out  using  a  grid  of  512  x 
200  points. 


Except  for  the  inclusion  of  pump  depletion,  the  present  version  of  BOUNCE  is  similar  to  the  one 
described  in  Ref.  (22).  (The  only  other  major  difference  is  that  the  predictor-corrector  algorithm  has 
been  replaced  by  a  split  operator  technique23  in  order  to  decrease  the  running  time.)  The  inclusion  of 
pump  depletion  is  a  nontrivial  problem,  however,  because  the  two  coupled  waves  are  incident  on  oppo¬ 
site  sides  of  the  Briilouin  medium,  and  neither  one  is  known  within  that  medium  at  the  outset.  One 
must  know  Es(x,z )  in  order  to  calculate  EL(x,z ),  and  vice  versa.  The  brute  force  solution  would  be  to 
simply  integrate  the  time  dependent  equations  over  many  optical  transits  across  the  medium  until 
steady  state  is  finally  attained.  The  approach  followed  here  is  a  modified  iteration  procedure,  which  is 
significantly  faster.  Basically,  it  starts  with  an  undepleted  pump  wave  E[l)  (x,z),  calculates  the  resultant 
backscatter  £$°  (x,z)  throughout  the  medium,  then  alternately  recalculates  these  waves  until  the  solu¬ 
tions  converge  to  self  consistency.  This  procedure  must  be  supplemented  by  stabilization  techniques  in 
order  to  achieve  that  convergence.  Acting  alone,  it  would  produce  solutions  that  oscillate  between  two 
extremes,  even  when  the  incident  pump  amplitude  is  "turned  on"  gradually  over  many  iterations. 


The  primary  stabilization  technique  is  based  on  the  fact  that  the  correct  steady  state  solutions 
must  satisfy  the  energy  conservation  relation. 


PL(z)  -  Ps(z)  -  PL(zO  -  />s(z,)  =  PL(zx), 


(2) 


where 

Pa(z)  -  S*'*2\Ea(x.z)\2  dz  (3) 

is  proportional  to  the  total  power  of  the  pump  (a  -  L)  or  backscatter  (a  «  5)  at  z.  Expression  (2)  fol¬ 
lows  immediately  from  (la,b)  and  the  stated  boundary  conditions  at  ±X/2.  The  actual  stabilizing  tech¬ 
nique  then  consists  of  multiplying  the  coupling  coefficient  g  in  Eqs.  (la,b)  by  the  factor 
1  +  yFiz)  -  y,  where 


F(z ) 


PL(z)  -  Ps(z) 
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PL(z)-Ps(z) 
PL(z{)  -  Ps(zJ 


(4) 


and  y  is  a  positive  bias  factor  of  order  3  <  y  <6.  Thus,  if  [^(jc.rj) I2  starts  to  become  too  large,  it 
automatically  limits  itself  by  reducing  the  growth  rate  1/2(1  +  yF  —  y)g\EL\2  in  Eq.  (lb).  A  similar 
restraint  prevents  excessive  depletion  of  the  counterpropagating  pump  wave  in  Eq.  (la).  When  the 
iterations  finally  converge  to  steady  state,  F  approaches  unity  throughout  the  medium,  so  the  bias 
becomes  inoperative.  In  a  one  dimensional  test  problem,  this  technique  worked  well  in  all  cases,  giving 
the  well  known  steady  state  solutions24  to  a  high  degree  of  accuracy.  In  two  dimensions,  however,  it 
failed  to  completely  stabilize  for  reflectivities  >  15%,  and  therefore  had  to  be  supplemented  by  an  addi¬ 
tional  technique. 


The  second  technique  performs  a  weighted  average  of  the  calculated  intensities  with  those  of  the 
previous  iteration.  In  the  «nh  iterative  solution  of  Eq.  (la),  one  first  calculates  \E[m]  (x,z)|2,  then 
makes  the  replacement 

\E[mHx,z)\2  —  1/2  \\E[m)  ix,z)\2  +  |£/m-1,(x,z) |2  P{Lm){z)l Pim~u(z)}  (5a) 

at  each  point  when  m  >  2.  A  similar  replacement  is  performed  when  solving  (lb): 

|£s,'">(x.z)l2  -  1/2  {|£sm,(x,z)l2  +  |£jm-1)(x,r)|2  Ps,w'(z)/^"-u(z)]  (5b) 

These  averages  have  the  important  property  that  they  leave  the  integrals  P[m'(z)  and  Psm>(: )  invari¬ 
ant,  and  therefore  do  not  interfere  with  the  primary  stabilization.  With  this  combination  of  techniques, 
<  the  solutions  typically  converge  to  self  consistency  within  15  to  20  iterations.  Best  results  were 

obtained  when  the  incident  £/""(x.z2)  was  turned  on  gradually  over  the  first  5  to  10  iterations. 
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In  addition  to  calculating  the  pump  and  backscatter  profiles,  BOUNCE  also  evaluates  the  average 
intensity  gain  6y  —  <|£s|2>z2/  <  |£s|2>:1,22  total  reflectivity  R  -  Ps(r2)/Pt(z2),  and  the  normalized 
correlation  function 

Hil)  ‘  P[G )fs<z)  \S-m  El (“)  EsUz)  *f '  (W 

which  measures  the  phase  conjugate  fidelity.  (Note  that  H  -  1  in  the  case  of  perfect  conjugation, 
£5  oc  El*.)  From  Eqs.  (la,b)  and  the  boundary  conditions  at  ±X/2,  one  can  readily  verify  that  H 
remains  constant  when  g  —  0;  thus,  H  *■  //(z2)  at  all  points  beyond  the  Brillouin  medium. 

III.  RESULTS 

Figure  1,  which  corresponds  to  the  case  of  negligible  pump  depletion,  models  a  typical  aberration 
compensation  experiment,  and  illustrates  several  general  features  of  the  simulations  presented  in  this 
paper.  A  coherent  hyper-Gaussian  pump  beam  of  X  -  694  nm  traverses  the  aberration  plate  A,  which 
imposes  a  random  phase  modulation  and  concomitant  —6  mrad  angular  divergence  ( — 25  x  diffraction 
limit)  shown  in  Fig.  2.  A  6  cm  lens  then  focuses  the  aberrated  beam  at  the  center  of  a  2  cm  thick  Bril¬ 
louin  medium.  In  all  of  these  simulations,  the  noise  source  at  was  modelled  by  a  random  complex 
field  Es(x,z |)  of  250  mrad  angular  divergence,  giving  typical  pump-noise  correlations  H(zx)  =  0.1%. 
However,  the  phase  conjugate  behavior  of  the  backscatter  was  found  to  be  generally  insensitive  to  the 
statistical  properties  of  the  noise  source  (including  even  a  coherent  source)  as  long  as  H(z{)  remained 
small. 


Referring  again  to  Fig.  1,  we  see  that  as  Es(x,z)  propagates  toward  z2,  it  develops  an  inhomogene¬ 
ous  intensity  distribution  similar  to  the  center  portion  of  |£i0c,z)|2.  This  tendency  of  the  pump  beam 
to  pull  the  backscattered  radiation  toward  the  strong  central  region  (i.e.,  spatial  gain  narrowing)  has 
been  noted  earlier,15- 16-'9- 22  and  appears  to  be  the  main  reason  for  the  lower  conjugate  fidelities  found 
in  focused  experiments.4*7  In  the  example  shown  here,  H  -  61%. 

From  z2,  the  backscatter  retraces  the  pump  path  through  the  lens  and  aberrator  to  produce  the 
near  and  far  field  intensity  profiles  shown  on  the  right  hand  side  of  Fig.  1.  The  phase  compensation  is 
evident  in  the  far  field,  where  the  backscatter  exhibits  a  strong  central  spike  that  closely  matches  the 
angular  spectrum  of  the  unaberrated  incident  beam.  The  —40%  unconjugated  component  appears 
mainly  as  low  level  sidelobes  and  background  hash;  hence,  expression  (6)  gives 

ff  ^  backscatter  power  in  the  central  spike  ^ 

total  backscatter  power 

From  a  rough  estimate  of  the  total  background  energy  shown  in  the  figure,  one  obtains  1  -  H  —  30%, 
with  the  remaining  10%  presumably  accounted  for  by  the  slight  broadening  of  the  central  spike,  plus 
residual  hash  lying  outside  the  ±4  mrad  interval  or  too  small  to  be  seen.  Although  the  peak  levels  of 
the  background  are  only  a  few  percent  of  the  on-axis  intensity,  they  are  still  significantly  larger  than  the 
experimentally  observed  values.4*7  This  stems  from  the  fact  that  the  hash  is  confined  to  only  one 
transverse  dimension;  i.e.  it  lacks  the  usual  geometrical  weighting  factor  that  would  require  much 
lower  levels  in  the  analogous  two  dimensional  "halo"  in  order  to  satisfy  expression  (7). 

In  the  near  field  of  the  aberrator,  the  conjugated  and  unconjugated  components  become  mixed, 
resulting  in  serious  degradation  of  the  entire  profile.  This  behavior,  which  was  found  in  all  cases,  has 
gone  largely  unnoticed  in  the  aberration  correction  experiments,  since  these  generally  record  only  far 
field  intensities.  It  will  be  discussed  in  greater  detail  in  Sec.  IV. 


Additional  calculations  were  carried  out  with  the  same  configuration  to  determine  how  the  perfor¬ 
mance  was  affected  by  the  gain  G /  in  the  absence  of  pump  depletion.  The  resulting  beam  profiles  and 


correlations  were  similar  to  those  shown  in  Fig.  1.  The  fidelity  exhibited  a  maximum  value  Hm3X  =* 
64%  around  G y  =s  e13,  and  decreased  gradually  to  55%  at  Gy  =  e25. 

The  elfect  of  pump  depletion  was  studied  by  increasing  the  average  noise  and  pump  intensities, 
while  holding  the  gain  approximately  constant  around  e17.*  Under  these  conditions,  the  fidelity 
improved  monotonicaily  with  reflectivity  R:  e..g.,  the  simulations  gave  H  «  68%  at  R  -  10%,  H  —  76% 
at  R  *  24%,  and  H  —  84%  at  R  —  47%.  Figure  3  shows  the  beam  profiles  for  the  last  case.  The  rea¬ 
son  for  the  improvement  is  evident  from  a  comparison  of  the  z2  plots  in  Figs.  1  and  3;  i.e.,  the  large 
backscatter  tends  to  counteract  spatial  gain  narrowing  by  selectively  depleting  the  more  intense  center 
portion  of  the  pump  beam  at  z  <  z2.  These  results  are  in  qualitative  agreement  with  a  recent  experi¬ 
mental  study  by  Mays  and  Lysiak,7  which  reported  a  steady  improvement  in  fidelity  with  increasing 
incident  power  well  above  threshold.  The  improvement  was  especially  evident  in  the  case  where  there 
were  no  reflections  from  the  walls  of  the  Briliouin  cell,  and  hence  where  gain  narrowing  effects  were 
most  important.  (In  two  other  experiments,6,9  the  fidelity  was  found  to  deteriorate  with  increasing 
incident  power;  however,  these  experiments  were  carried  out  under  conditions  where  self  focusing  is 
likely  to  have  dominated  within  the  cell.) 

Although  the  enhanced  fidelity  significantly  reduces  the  background  hash,  and  narrows  the  axial 
spike  at  the  far  field  (as  one  would  expect  from  expression  (7)),  it  results  in  only  a  partial  improvement 
in  the  near  field  beam  quality  at  the  aberrator.  Figure  4  shows  that  in  spite  of  the  persistent  small  scale 
inhomogeneities,  the  backscatter  can  still  retrace  gross  features  imposed  upon  the  incident  beam.  This 
result  is  consistent  with  the  image  replication  experiment  of  Ref.  (11),  in  which  an  aberrator  was 
present  in  addition  to  the  mask.  It  is  interesting  to  note  from  the  z2  plot  that  the  effect  does  not 
require  an  image  of  the  structure  to  be  maintained  within  the  Briliouin  medium. 

In  all  of  the  simulations  presented  here,  the  aberrator  was  modelled  using  the  same  statistical  real¬ 
ization.  Other  realizations  yielded  different  values  of  H  and  differences  in  detailed  structure  of  the 
profiles;  however,  the  trends  and  general  appearance  of  the  profiles  remained  unchanged.  In  particular, 
the  poor  beam  quality  at  the  aberrator  was  found  in  all  cases,  even  one  for  which  the  fidelity  was  88%. 

IV.  DISCUSSION 

Several  of  the  results  presented  in  Sec.  Ill  were  not  unexpected;  e.g.,  the  enhancement  of  fidelity 
by  the  pump  depletion  was  suggested  earlier  by  qualitative  theoretical  arguments.15  19,22  However,  the 
consistently  poor  beam  quality  of  the  backscatter  at  the  aberrator,  even  for  H  =  90%,  deserves  some 
additional  comment.  This  behavior  can  be  adequately  explained  if  one  approximates  the  pump  and 
backscatter  amplitudes  at  the  right  hand  side  of  the  aberrator  by 

El  =  constant  (8a) 

£,(*)  =  Rx,2El\  1  +  £  exp  (2ir/A&f/D)]  (8b) 

l  <v*o  J 

within  the  interval  -D/2  <  x  <  D/2,  and  zero  outside.  The  mode  amplitude  can  then  be  written  as 
ay  —  Esn/ESq,  where  £so  and  Ess  are.  respectively,  the  peak  far  field  amplitudes  of  the  central  spike 
and  the  Nth  lobe  at  angle  Nk/D.  Substituting  Eqs.  (8a, b)  into  (6),  one  obtains 

1  +  £  I  at  v  1 2  (9) 

i»j*0 

which  is  consistent  with  expression  (7)  if  all  lobes  are  of  equal  width.  For  M  modes  of  equal  RMS 
amplitude  a*  1/5,  the  fidelity  is  H  =  (l  +  A/aaMS)-1  while  the  fractional  fluctuation  of  |£s(x)|2  is  of 
order  ±2A/1/2  a*,ws  (assuming  2A/,/2  arRMS  <  1  and  aRMs  «  D-  For  the  simulation  shown  in  Fig.  3, 
we  estimate  M  =s  10  and  a^Ms  28  0.01.  obtaining  H  =  91%  with  4:1  intensity  variations  at  the  aberra¬ 
tor.  It  is  instructive  to  note  that  even  a  99%  fidelity  would  result  in  ±20%  variations. 
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The  simulations  presented  in  this  paper  have  shown  the  importance  of  pump  depletion  and  spatial 
gain  narrowing  on  phase  conjugation  in  focused  SBS  experiments,  and  illustrated  the  problem  of  near 
field  beam  quality.  A  modified  version  of  the  numerical  techniques  described  here  can  also  be  used  to 
model  an  optical  waveguide  Brillouin  mirror,  or  a  degenerate  four  wave  mixing  geometry  in  the  limit  of 
small  beam  crossing  angles. 
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Fig.  1  —  SBS  phase  conjugation  model,  showing  intensity  profiles  of  the  incident  light  (solid  lines) 
and  backscatter  (dotted  lines)  for  the  case  of  negligible  pump  depletion 
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